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Conservation laws for the Maxwell-Dirac equations with a dual Ohm's law 

Nail H. IbragimovQand Raisa Khamitova 

Department of Mathematics and Science, Research Centre ALGA: Advances in Lie Group Analysis, Blekinge Institute 
of Technology, SE-371 79 Karlskrona, Sweden 

Bo Thidffl 

Swedish Institute of Space Physics, Angstrom Laboratory, P. O. Box 537, SE-751 21 Uppsala, Sweden^ 



o 
o 



in 



> ■ 
oo ■ 

o ■ 
m ■ 
O ' 
r» . 
o ■ 

^ : 

Of 



Using a general theorem on conservation laws for arbitrary differential equations proved by Ibragimov, we have 
derived conservation laws for Dirac's symmetrized Maxwell-Lorentz equations under the assumption that both 
the electric and magnetic charges obey linear conductivity laws (dual Ohm's law). We find that this linear system 
allows for conservation laws which are non-local in time. 



PACS numbers: 1 1 .3O.-i,O2.20.-a,02.3O.Jr,03.5O.De,41 .20.Jb 



I. INTRODUCTION 

In all areas of physics, conservation laws are essential since 
they allow us to draw conclusions of a physical system under 
study in an efficient way. 

Electrodynamics, in terms of the standard Maxwell electro- 
magnetic equations for fields in vacuum, exhibit a rich set of 
symmetries to which conserved quantities are associated. Re- 
cently, there has been a renewed interes t in the utilisation of 
such quantities. Here we use a theorem of Ibragimov (2006) to 
derive conservation laws for Dirac's symmetric version of the 
Maxwell-Lorentz microscopic equations, allowing for mag- 
netic charges and magnetic currents, where the latter, just as 
electric currents, are assumed to be described by a linear rela- 
tionship betwee n the field and the current, i.e., an Ohm's law. 
The method of llbragimovl d2006l) produces two new adjoint 
vector fields which fulfil Maxwell-like equations. In particu- 
lar, we obtain conservation laws for the symmetrized electro- 
magnetic field which are non-local in time. 



We employ the usual convention of summation in repeated 
indices. 

Recall that a necessary condition for extrema of a varia- 
tional integral 



£j(x, M,M(i))dx 



(1) 



with a Lagrangian L(x,u,u^), depending on first-order 
derivatives, is given by the Euler-Lagrange equations 

8L dL ( dL\ 

^^- D <W H a=1 '-' m - (2) 



We will understand by a symmetry of a certain system of 
differential equations a generator 



X = % i (x,u) 1 r-; + r]< x (x,u) 1 — 
ox' our 



(3) 



X 



II. PRELIMINARIES 
A. Notation 

We wil l use the following notation (see, e.g., Ref. 
Ibragimov, Il999f) . Let x — (x 1 , . . . , x") be independent vari- 
ables and u — (u l , . . . ,u' n ) be dependent variables. The set 
of the first-order partial derivatives uf — du a /dx l will be 

denoted by ui\\ — {uf}, where a = 1, ,m and = 

1, . . . ,n. The symbol Dj denotes the total differentiation with 
respect to the variable x l : 



D; 



d_ 

"ox 1 



du a 



11 duf 



u a 3 
ijk du% 
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of a continuous transformation group admitted by differential 
equations under consideration. 

A vector field C = (C 1 ,. . . ,C") is said to be a conserved 
vector for the differential equations (fJJ if the equation 



Di a =o 



(4) 



holds for any solution of Eq. (0. 

If one of the independent variables is time, e.g., x" — t, then 
the conservation law is often written in the form 



dE 
dT 



= 0, 



where 



E= J C n (x,u(x),u {1) (x))dx 1 ■■■dx"- 1 . (5) 
Accordingly, C" is termed the density of the conservation law. 
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B. Basic conservation theorem 



We will employ the recent general theorem (Ibragimov, 
2006) on a connection between symmetries and conserva- 
tion laws for arbitrary systems of sth-order partial differential 
equations 



then Eqs. (O provide r conserved vectors C\ , 
components 



,C r with the 



M = l, 



,r; i 



F a (x,u,U( 



>"(i))=0, a = \. 



,m, 



(6) 



ELECTROMAGNETIC EQUATIONS 



where F a (x,M,M (1) , . . . ,u {s) ) involves n independent variables A Basic equations and the Lagrangian 

x = (x l , . . . ,x n ) and m dependent variables u — (u l , . . . , u m ), 



u = u(x) together with their derivatives up to an arbitrary order 
s. For our purposes, we formulate the theorem in the case of 
systems of first-order dif ferential equation s. 
Theorem ELI. (See Ref. llbragimovll2006l Theorem 3.5). Let 
an operator (O be a symmetry of a system of first-order partial 
differential equations 



F a (x,u,uii) 



0, a=\,. 



where v — (v 1 , . . . , v m ). Then the quantities 



C 



i=\. 



(7) 



(8) 



furnish a conserved vector C = (C 1 , . . . ,C n ) for the equations 
d2]i considered together with the adjoint system 



<5£ 



F*(x,u,v,u w ,v (l) ) = — =0, a = l 



(9) 



where 



8u c 

8 d d 
8u a ~ du a ~ D< du?> " ' 



and v = (v\...,v m ) are are new dependent variables, i.e., 

v = v(x). 

Remark DLL The simultaneous system of equations (0 
and (O with 2m dependent variables u = (u , . . . ,u m ), 



'=(v 



. ,v m ) can be obtained as the Euler-Lagrange equa- 



tions (O with the Lagrangian 

£ = ^Fp(x, U, «(,)) 



Indeed, 



= F a (x,u,ur 



5£ 
<5iL 

^ 5 =F a (jf,II,V,«(i),V 



(i) 



(10) 

(ID 
(12) 



Remark II.2. The conserved quantities ([H) can be written in 
terms of the Lagrangian (TTOb as follows: 



Remark II.3. If Eqs. (0 have r symmetries X\ , . . 
form (01, 

X fl = ^ fl (x, U )^+r]«(x )U )^, fl = l, 



(13) 
,X r of the 

,.,r, 



Adopting Dirac's ide as on the existence of magnetic 
monopoles dPirad, 1 193 lb . one can formulate a s ymmetrized 
versio n of Maxwell's electromagnetic equations dSchwingen 
119691) . In SI units an d in microsco pic (Lorentz) form, these 
equations are [cf. Ref. Thide, 2006, Eqs. (1.50)]: 



VxE + — + jUoj,„ = 0, 
1 <?E 

VxB--^ — - ju j e = < 

V-E-^ c 2 p e = 0, 
VB - Hop m = 0, 

together with the dual Ohm's law 

Se = C7 e E, j,„ = <T,„B. 



(14a) 

(14b) 

(14c) 
(14d) 



(15) 



where a m and a e are constant scalar (rank zero) quanti- 
ties. The first equation in ( TT3T > is Ohm's law for electric 
currents. The second equation is a dual Ohm's law for 
magnetic curr ents, that, fo r symmetry reasons, was intro- 
duc ed in Ref. Thide, 200g, Eqs. (2.60); see also Eq. (5 ) in 
Ref. iMever- Vernetl 1 19821 Eg. ( 38) in Ref.l01esenlll996l and 
its generalization Eq. (8) in Ref. lCoceal et all 1 19961 

Now we substitute Eqs. ( fl3T > into Eqs. d!4a| )-( fT4bl i. The en- 
suing equations involve, along with the light velocity c, three 
other constants, a e , G m and /io- We eliminate two constants by 
setting 



~ Mo 
t = ct, B = cB, a e = c[ioo e , a m = — a m , 

c 



(16) 



and rewrite our basic Maxwell-Dirac equations ( TBI , discard- 
ing tilde, as follows: 



c)\\ 

VxE + — - 

dt 



VxB 



a m B 

0,E = 



= 0, 
0. 



(17) 



dE 

~dt 

V-E-p e = 0, 
VB-p m = 0, 

The system ([T7| > has eight equations for eight dependent vari- 
ables: six coordinates of the electric and magnetic vector 
fields E = (E 1 ,E 2 ,E 3 ) and B = (B 1 ,B 2 ,B 3 ), respectively, and 
two scalar quantities, viz., the electric and magnetic monopole 
charge densities p e and p,„. 
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Using the method of llbragimovl d2006l) we write the La- 
grangian ( fTQb for Eqs. ( fT7] > in the following form: 



dB 

L =V- (VxE+-^- + a m B) +^(V-E-p, 

W (VxB-^-o,E) + /?„,( V-B-p,, 



(18) 



where V,W,/? e ,7? m are adjoint variables (we note in passing 
that V is a pseudovector and R,„ a pseudoscalar). With this 
Lagrangian we have: 



<5£ 
8V 


= VxE + 


dB 




<5£ 


<5£ 
<5W 


= V xB 


<?E 

~ d7 


-cr e E, 


<5iL 
5K« 



(19) 



V-B-p m , 



and 



— =VxV+^-o e W-Vfi e , — = -R e 
oE o*f dp e 



— = V x W — + a m \ - VR, 



8L 



(20) 



<5B dt ' 8p„ 



R„ 



It follows from Eqs. (fT9b — (f20b that the Euler-Lagrange equa- 
tions (12) for the Lagrangian ( TT8l provide the electromagnetic 
equations ( fTTI i and the following adjoint equations for the new 
dependent variables V, W,R e ,R m : 



dW 

Vx V + — -a e W = 0, 
dt 

dy 

V x W - — + a m V = 0, 
dt 

Re = 0, R m = 0. 



(21) 



Remark III.l. Let the spatial coordinates x l ,x 2 ,x 3 be x,y,z. 
For computing the variational derivatives SH / 8E and SL/SB 
in Eqs. ( f20l >. it is convenient to use the coordinate representa- 
tion of the Lagrangian ( fT8l . namely: 

£ = V 1 (£ y 3 - £? + B] + a m B x ) + V 2 (E\ - E 3 + B 2 + o m B 2 ) 
+ V 3 [E 2 -E]+ B 3 + a, n B 3 ) + R e {E\ + E 2 + E 3 - p e ) 

(22) 

+ W 1 (By —B\ — E] - a e E l ) + W 2 (B\ -B 3 x -E 2 - a e E 2 ) 
+ W 3 (B 2 - B l y - E 3 - o e E 3 ) + R m {B\ + Bl + B 3 -p m ). 



B. Symmetries 

Eqs. ( fTTb are invariant under the translations of time t and 
the position vector x = (x,y,z) as well as the simultaneous 
rotations of the vectors x,E and B due to the vector formu- 
lation of Eqs. ( fTTI i. These geometric transformations provide 



the following seven infinitesimal symmetries: 





d 

X = Ft , X, 


d 

dx 1 


x 2 -- 


d 


d 

X3 = Tz> 






d d 

^ dx dy 


- F z 


-E 


d 

dE 2 ' 


VB 2 — r —B 

dB 1 


d 

dB 2 ' 




d d 

dx dz 


, d 

-E 3 ^—r 

dE l 


-E l 


d 

dEl^ 


-B 3 ^—r-B l 
dB 1 


d 


X23 


d d 

Z dy ^ dz 


, d 
dE 2 


-E 2 


d 

dE 1 ^ 


id 

-B 3 —r-B 2 

dB 2 


d 

dB 1 ' 



(23) 

The infinitesimal symmetries for the adjoint system (BTI) are 
obtained from (|23l by replacing the vectors E and B by V and 
W, respectively. Moreover, since Eqs. ( fTTI i are homogeneous, 
they admit simultaneous dilations of all dependent variables 
with the generator 



: E • 



dE 



B 



dB 



Pe 



dpe 



dp,, 



(24) 



where 



d 3 j d d 3 ,. d 

E -dE-h 1 E dE^' dB't*™' 



dB 1 



Recall that the Maxwell equations in vacuum admit also 
the one-parameter group of Heaviside-Larmor-Rainich dual- 
ity transformations 



E = Ecosa — Bsina, 
with the generator 



B = Esina +Bcosa 



(25) 



X 



B'- 



dB i dE* 

Also recall that the "mixing angle" a in (|25| ) is a pseudoscalar. 

It was shown in jlbragimovi 12006) that the group ( f25l l pro- 
vides the conservation of energy for the Maxwell equations. 
Let us clarify whether Eqs. dTTb admit a similar group. Let 
therefore 

d d d d 
X = E- — B — +p e - p m3 — . (26) 

dp m dp e 

1 is written 



<9B <5E 
The prolongation of the operator 

d 



X=E- 



dB 

E,- 

-Ey 



B 


d 


d d 






dE + 


dp m dpe 






d 


B, 


d d 

hE r 

dE, x dB x 




d 


dB, 


d~E~ x 


d 

"d~B^, 


-By 


dEy Z dB, 


B z 


d 

"d~E z 



(27) 



Reckoning shows that the operator (f27j acts on the left-hand 
sides of Eqs. ( flTT i as follows: 

X ( V x E + B, + a m B) = - ( V x B - E, - a, n E) , 

l(VxB-E,-ff f E) = VxE+B f + a f B, 

X(V-E-p e ) =-(V.B-p m ), 

X(V-B-pb) = V-E-p e . 
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It follows that the operator d26l i is admitted by Eqs. ( TP71 ) only 
in the case 



a„, = a,. 



(28) 



IV. CONSERVATION LAWS 

A. Derivation of conservation laws 

We will write the conservation law in the form 

A(T)+div* = 0, (29) 

where the pseudoscalar T is the density of the conservation 
law d29l ), the pseudovector current % = (x 1 ,X 2 >X 3 )> and 

div* = V x = D x (x l )+D y {x 2 )+D z {x 3 ). 

Let us find the conservation law furnished by the symme- 
try ( |26] i when the condition d28| i is satisfied, a,„ = O e . Ap- 
plying the formula <[8J to the symmetry (|2oT i and to the La- 
grangian (fT8l . we obtain the following density of the conser- 
vation law i 



Thus, 



dL dH 

T = E • — B — = E V + B W. 

<9B, <?E, 



T = E-V + B-W. 



(30) 



The pseudovector % is obtained likewise. For example, using 
the Lagrangian in the form d22i >. we have: 



A dB z dE x 



The other coordinates of % are computed likewise, and the 
final result is: 



£ = (ExW)-(BxV). 



(31) 



One can readily verify that (f30b and OTb provide a conser- 
vation law for Eqs. ( fTvT > considered together with the ad- 
joint equations (f2TT >. Indeed, using the well-known formula 
V ■ (a x b) = b • (V x a) - a • (V x b) and Eqs. <fJ3 and (ED, 
we obtain: 

A(t) =E ( -V + E-V, + B,-W + B-W, 

= V (V x B a e E) +E- (V x W + a m V) 
— W • (V x E + c m B) — B • (V x V — oVW), 

V-£=V-(ExW)-V-BxV) 
= W • (V x E) - E • (V x W) 
-V-(VxB)+B-(Vx V). 

Whence, 

D, (t) + div* = (cj,„ - a e )(E • V - B • W). 

It follows again that the conservation law is valid only if 
a m = o e . 



Remark IV.l. The conservation law given by (T3Qb— OTb de- 
pends on solutions (V,W) of the adjoint system ( f2Tb . How- 
ever, substituting into Eqs. ( 130b and OTb any particular solu- 
tion (V, W) of the adjoint system OTb with a m = a e , one ob- 
tains the conservation law for Eqs. ( fTTb not involving V and 
W. Let us denote a m = a e = a and take, e.g., the following 
simple solution of the adjoint system (f2TJ): 



V l =, 



V 1 



y j = 0; w l 



w = w 3 



0. 



Then Eqs. OTJ-tEB yield: 



T=(E 1 +B l )e at 
X '=0 X 2 = (E 3 -B 3 )e a > X 3 = (B 2 -E 2 )e at . 



Remark IV.2. The operator (1261 ) generates the one-parameter 
group 

E = Ecosa — Bsina, B = Esina+Bcosa, 
p e = p e cosa — p,„sina, p '" = p e sin a + p m cos a . 

where, again, the "mixing angle" a is a pseudoscalar. 
Remark IV.3. In the original variables used in Eqs. J 14-ab — 
(1 14bb and dT5l >. the operator d26l i is written: 

Id d did 

X = -E-^-cB- — +cp e - Pm-5— • 

c dB dE dp m c dp e 

Applying similar calculations to the generator ( |24| ) of the 
dilation group provides the conservation law with 

t = B V E W, £ = (ExV) + (BxW). (32) 

This conservation law is valid for arbitrary a,„ and o e . Indeed, 

D t (z) = B,-V + B-V,-E,-W-E-W f 

= -V-(VxE + a m B)+B-(Vx W + a„,V) 

- W • (V x B — a e E) + E • (V x V - c e W) 
= -V-(VxE)+B-(VxW) 

-W-(VxB)+E-(Vx V), 

V • % = V • (E x V) + V • (B x W) 
= V-(VxE)-E-(Vx V) 
+ W-(VxB)-B-(VxW). 

Hence, D, (x) + div£ = 0. 

Let us find the conservation law provided by the symmetry 
X Q = d jdt from (|23). Formula © yields: 

T = £-E~-B t ~=£+E r W-B r V. 

oE t dB t 

Since the Lagrangian L given by ( fT8l vanishes on the solu- 
tions of Eqs. (fTTT i, we can take 



T = E,-W-B,-V 



(33) 



or 



t = W • [(V x B) - a e E) + V • [(V x E) + o m B]. (34) 
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Let us calculate the pseudovector %. Formula © yields: 

X = E ' dE x ' d^~ x 
Using the Lagrangian in the form d22i l. we have: 

X 1 = -E 2 V 3 + E 3 V 2 - Bjw 3 + B 3 W 2 . 

The other coordinates of % are computed similarly, and the 
final result is 



X = (V x E,) + (W x B, 



(35) 



Thus, the time translational invariance of Eqs. ( fT7b leads to 
the conservation law d29j > with T and X given by ( [34| > and fl35l >, 
respectively. 

Remark IV.4. Let us substitute in Eqs. (l34t and d35l l the fol- 
lowing simple solution of the adjoint system (cf. Remark 
BED: 

yl = e Omt^ y 2 = y3 = q. ^1 = g (T e f ^ ^2 = ^3 = q 

Then Eqs. d34t-(l35l) yield: 



T = (fiy - B\ - a e E 1 + (E$ - + a m B 1 )e (7 '" r , 
Z 1 = 0, X 2 = -E?e a "'' - B 3 e a ?' , ^ 3 = E 2 e a '"' + B 2 e a '' . 

The conservation law provided by the symmetry X\=d / dx 
from ( |23l has the following density: 

dL <?£ 
T = -E.v • — - B, • — = E. r • W - B x • V . 

<7E, dB t 

For the pseudovector £ the formula ^ yields: 
i dL d£j 

Using the Lagrangian in the form d22l l. we have: 

X 1 = £ - £ 2 V 3 + £ 3 V 2 - # 2 W 3 + B 3 W 2 . 

The other coordinates of X are calculated similarly: 

X 2 = E\v 3 -E 3 X V X +B X W 3 -B 3 .W 1 , 
X 3 = -E^+E^V 1 -B x W 2 +B 2 x W l . 

We can ignore L in X 1 since = on solutions of Eqs. (ITTb 
and (f2TT >. and the final result is 



r = (V X E,) + (W x B,) 



(36) 



Replacing x by y and z we obtain the following conservation 
laws corresponding to X2 — d/dy and X3 = d/dz, respec- 
tively: 



and 



T = E y • W — By -V, X = (V x E v ) + (W x B Y ) 



T = E Z -W-B Z -V, Z = (VxE z ) + (Wx B, 



Applying Formula (0 to the symmetry X\2 and to the La- 
grangian (l22l >. we obtain the following density of the conser- 
vation law: 



=W 2 E 1 - W X E 2 + (yE x -xEy) ■ W 
- (y 2 ^ 1 - V l B 2 ) - (yB x -xB y ) • V. 
The densities of the conservation laws forXi3 and X23 are 

t = W 3 E l - W l E 3 + (zE x -xE z ) ■ W 



- {V 3 B l - V - (zB x -xB z ) • V 



,1 D 3 



and 



X = W 3 E 2 - W 2 E 3 + (zE v - yE z ) ■ W 

- (V 3 B 2 - V 2 B 3 ) - ( Z B y -yB z ) • V, 

respectively. Finally, the densities of conservation laws corre- 
sponding to the rotation generators Xjj can be written as one 
vector: 

T = WxE + W-(xx V)E -VxB-V-(xx V)B, (37) 

where x = (x,y,z). 

The operator X12 provides the following pseudovector £: 

X 1 = -V 3 E l -y(E 2 V 3 -E 3 V 2 ) + x(E 2 V 3 -E 3 V 2 ) 



W 3 B l -y{B z W 3 -B S X W Z ) +x(B^W 3 -B s y W l ), 



32w3 d3ti/2 



2ii/3 d3ti/2n 



V 3 E 2 +y(E l x y i -E^V 1 ) -x(E x l V i -E*V 



,U,3 



,l T/ 3 



W 3 B 2 + y(B x W 3 - B X W 1 ) - x{B l W 3 - B 3 W 1 ) , 



>'w3 d3ti/1 



hi/3 d3ii/1x 



^ J =V 1 £ 1 +V Z E Z -y(E x V z -E z V l )+x(E y l V z -E z V l ) 

+ W l B l + W 2 B 2 -y{B l x W 2 - B 2 x W l ) +x(B^W 2 - B 2 ,W l ) . 

The pseudovector X f° r the operator X^ has the following 
form: 

X 1 =V 2 E l ~z{E 2 V 3 ~E 3 V 2 ) +x(E 2 V 3 -E 3 V 2 ) 

+ W 2 B l -z{B 2 x W 3 -B 3 x W 2 )+x(B 2 W 3 -B 3 W 2 ), 
X 2 = -V l E l - V 3 E 3 + z{E\ V 3 - E 3 V 1 ) - x{E\V 3 - E 3 V 1 ) 

- W x B l - W 3 B 3 +z{B l x W 3 ~B 3 x W l ) -x{B\W 3 -B 3 z W l ), 
X 3 =V 2 E 3 - z{E l x V 2 -E 2 V l )+x(EW 2 -E 2 V v ) 

+ W 2 B 3 -z{B x W 2 - B 2 x W 1 ) + x(B] W 2 - B 2 W 1 ) . 
The operator X23 provides the following pseudovector X- 
X 1 =V 2 E 2 + V 3 E 3 -z{E 2 V 3 -EyV 2 ) +y{E 2 V 3 -E 3 Z V 2 ) 

+ W 2 B 2 + W 3 B 3 - z{B 2 W 3 - B 3 W 2 ) +y{B 2 W 3 - B 3 Z W 2 ), 
X 2 =-V 1 E 2 +z(E y l V 3 -E^V 1 ) -y{E}V 3 -E 3 V l ) 

-W x B 2 + z{b]MV 3 -B 3 ,W l )-y(B l z W 3 -B]w x ), 

x 3 



" ; - -V l E 3 -z(E y l V 2 -E 2 V l )+y(ElV 1 -E 2 V 1 ) 
- W l B 3 - z{B y W 2 - B 2 W 1 ) +y{Bl W 2 - B 2 W l ) . 
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B. Two-solution representation of conservation laws 



The conserved quantities obtained in Subsection IIV.AI in- 
volve solutions V, W of the adjoint equations (fJTJ. It may be 
useful for applications to give an alternative representation of 
the conserved quantities in terms of the electric and magnetic 
vector fields E,B only. 

We suggest here one possibility based on the observation 
that one can satisfy the adjoint system (fJTJ by letting 



V(v)=B(x,-f), 

W(v) = E(x,-f), 

R e (x,t) = V-E(x,-f)-p e (x,-f), 

R m (x,t) = V-B(x,-f)-p m (x,-f), 



(38) 



where E(x,i),B(x,i) solve Eqs. ([T7l > with s = —t. Indeed, 
employing the substitution d38l l and the notation s = —t we 
have 

<?W 

Vx V+ — -cr e W 
at 

ds at 



Vx W--^ + cr m V 



(39) 



= V x E x,j t^-t + o m B{x,s) , 

as at 

R e = V-E(x,s)-p e (x,s) = 0, 
R m = V -B(x,s) - p m (x,s) . 

Hence, the adjoint equations OTT i reduce to ( fTvT >: 

_ (9B(x,i) „. 

VxEx,s+ — ^— ^ + a,„B (x,5 =0, 

(75 

V-B(x,s) ^-a e E(x,i)=0, (40) 

V-E(x,i)-p e (x,i) = 0, 
V-B(x,i)-p m (x,i)=0. 

Let (E(x,f),B(x,f)) and (E'(x,f ),B'(x,f )) be any two so- 
lutions of the electromagnetic equations ( fTTI ). Substituting in 



(El the solution (E',B'), we obtain the two-solution repre- 
sentations of the conservation laws. For example, the conser- 
vation law given by (f30b — (f3TT> has in this representation the 
following coordinates: 



T = E(x,f)-B'(x,-f)+B(x,f)-E'(x,-f), 

X = [E(x,r) x E'(x, -t )] - [B(x,f ) x B'(x, -t)} . 



(41) 



In particular, if the solutions (E(x,f ),B,f )) are identical, (I4T1) 
provides the one-solution representation: 



T = E(x,f ) ■ B(x, -t ) +B(x,f ) • E(x, -t), 

X = [E(x,f ) x E(x, -t)\ - [B(x,f ) x B(x, -t )] . 



(42) 



All other conservation laws can be treated likewise, e.g., 
the conservation law given by d33l and ( f35l > has the following 
two-solution representation: 



T = E,(x,r) -E'(s,-f) -B f (x,f) .B'(x,-f), 
X = [B'(x, -0 x E,(x,r)] + [E'(x, -t) x B,(x,f) 
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